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Abstract 

Entanglement in random states has turned into a useful approach to quantum 
thermalization and black hole physics. In this article, we refine and extend the ‘ran¬ 
dom unitaries framework’ to quantum field theories (QFT), and to include conserved 
charges. We show that in QFT, the connection between typical states, reduced sub¬ 
systems and thermal dynamics is more transparent within the Fock basis. We provide 
generic formulae for the typical reduced density matrices and entanglement entropies 
of any given subset of particles. To illustrate our methods, we apply the generic 
framework to the simplest but non trivial cases, a massless scalar field in two dimen¬ 
sions and its generalization to the case of N scalar fields, including the large N limit. 
We find the effective temperature, by matching the reduced dynamics to a Gibbs 
ensemble, and derive the equation of state of the QFT. The deviations from perfect 
thermality are shown to be of order 1/5 instead of exp(—5), a result which might 
be relevant for black hole physics. Finally we describe the analogue of the so-called 
‘Page curve’ in the QFT scenario as a function of the energy scale which divides high 
from low energy degrees of freedom. 
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1 Introduction and summary of results 


1.1 Quantum thermalization and random matrices 

In its simplest formulation, the black hole information paradox [T] follows from the state¬ 
ment that black hole radiation [2] is precisely thermal, with the temperature depending 
just on the conserved charges defining the black hole. Given that there is an exponentially 
big number = e Sw of initial states that could collapse to the given black hole, where 
S'bh is the Bekenstein-Hawking entropy M, it is concluded that after the black hole has 
evaporated the information about the initial state is lost. Mathematically, the paradox 
arises since under unitary evolution of the initial quantum state, we can never produce a 
thermal ensemble, 

f/|^)(^|G t ^p Gibbs , (1.1) 

where p Glbbs is an appropriate Gibbs distribution for the problem at hand. Interestingly, 
this problem is strikingly similar to one old problem in quantum mechanics, the problem 
of ‘quantum thermalization’, see SUSIE] and references therein. The problem of quantum 
thermalization can be easily stated as: 


In the context of microscopic unitarity, and given (1.1), how can Gibbs ensembles 


emerge: 


? 


In what follows we will assume the Hilbert space to be factorizablcQ : 




( 1 . 2 ) 


The number of factors in (1.2) is unconstrained, it can be finite or infinite. Besides, the 
Hilbert space factors do not need to be equal to each other. If the algebra of operators 
acting on Hi is given by O "', with a; a discrete parameter, the algebra of operators acting 
on H is 


o = o j 11 <g> <g> ■ ■ • o: 


(1,3) 


Reduced subalgebras/subsystems A are defined by specifying a set of Hilbert factors, and 
considering non-unity operators only in those subfactors: 


O a = O a <g> 1 <8> 1 <g> • • • . (1.4) 

As is well known, given a unitarily evolving pure state p{t) = |'0(f))('0(f)|, the expectation 
value of any operator belonging to such subalgebra A is completely characterized by its 

lr The ideas and results concerning the problem quantum thermalzation are most probably not sensitive 
to this condition. But in this framework the concepts can be explained more transparently, and a plethora 
of results can be derived. Besides, the models we will explicitly consider in this article fall in this factorizable 
class. 
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associated reduced density matrix: 


{ip(t)\0 A \ip(t)) = TV(p*(«)O a ) , 


(1.5) 


where 

PA(t) = Tr A (p(t)) , (1-6) 

and A is the complement subalgebra/subsystem of A. 

Although perfect thermality cannot be attained within unitary evolution, we might 
still expect certain approximate thermality for the correlation functions of the theory. The 
condition stating that all measurements (1.5) are approximately thermal 


(m\0 A \m) = Tr(p G ‘ bb, £> A ) ± error 


(1.7) 


can be more rigorously stated as 


PaW = P A bta ± error , 


( 1 . 8 ) 


with small errors, and where 

PA bb ‘ = Tr A (p Glbl ”) , (1.9) 

is the reduced thermal density matrix associated to A. 

Given these preliminaries, some immediate questions concerning quantum thermaliza- 
tion are: 


• What is the corresponding p Glbbs for the specific theory and state p(t) considered? 
Or equivalently, how the parameters defining p Glbbs are related to the parameters 
defining p{t)l 

• How does the error scale in the thermodynamic limit? 


Considering (1.8) and the previous questions, one modern approach to tackle this problem 
is to compute the entanglement entropy of the reduced subsystem A, which is defined by 


si(t) = -TV(pa(() log-p A (f)) . 


( 1 . 10 ) 


If our state is such that relation (1.8) holds, the next insightful relation follows as well: 

S£(t) = S'cibbs ± error , (1.11) 


which states that the origin of thermal entropy in quantum mechanics lies in quantum 
entanglement. 
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The discussion has been until now entirely generic. To proceed we should take a specific 
Hamiltonian H and find the exact evolved quantum state: 

p(t) = \m)W )I = tW W = e- m \i,)^\e- m . (1.12) 

This is certainly a difficult task, since we are usually interested in high energy states, with 
many particles excitecQ But due to these complicated dynamics, we expect that some 
approximations can be done. In particular, one common procedure is to assume that after 
a certain time we can approximate the evolved state as a random vector^] 

k t re ) ax )) — It/^random) i ( 1 . 13 ) 

which can be equivalently thought as being created by the application of a random unitary 
matrix to a given initial state: 


(^random) — ^random |Vv ■ 


( 1 . 14 ) 


The question now is obvious: 

• How do we define ^random) or U vandom 7 

If we could answer this question generically for any theory and initial state, it would allows 
us to compute typical values of physical quantities in a unitary framework. The physical 
quantities in which we are interested in this article are the typical reduced density matrix 
of a given subsystem A, given by 


pT* = JdpTrM , ( 1 . 15 ) 

and the typical entanglement entropies 

= JdpSi(p). ( 1 . 16 ) 

This approach was pioneered in Haig. In those articles a bipartite Hilbert space was 
considered, with dimension AB } where A and B are the dimensions of the corresponding 
factors of the total Hilbert space. For this Hilbert space, (^random) was defined by averaging 
over the continuum manifold of quantum states PI. with a constant measure, i.e the 

2 Some important progress has been done in the context of integrable systems, see [7| and references 
therein, conformal field theories, see 0, and in the context of holographic thermalization, see PHU]. In 
this article we will be interested in other approach, which is more general in certain aspects and it is able 
to provide the scaling of the deviations from thermality with the size of the system. 

3 See m and |I2| for two nice expanded reviews on the subject. 

4 We use the words ’’typical”and ’’average” interchangeably. For example, the typical reduced density 
is by definition the same as the average density matrix. 
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Haar ensemble of quantum states in the full Hilbert space. Again, this is equivalent to 
applying a random unitary U T an dom picked from the Haar ensemble of random unitaries. 
With this specific definition of randomness in the given Hilbert space, the typical reduced 
density matrix (1.15) was found in [X3J and, using those results, the typical entanglement 
entropy (1.16) was found in j5rl In the light of the previous questions about quantum 


thermalization, their results for this specific case are the following: 


• The typical reduced density matrix of a subsystem A is proportional to the identity 
matrix. Therefore, this reduced density matrix might be associated to the thermal 
density matrix of any Hamiltonian (defined for the finite dimensional Hilbert space) 
in the infinite temperature limit. 


The error depends on the size of both subsystems A and B. For A ~ 0( 1), the error 
in formula (1.11) is of O(j^) ~ 0(2~ s ), where S is the total entropy of the system. 
For A = B the error is of 0(1). 


The first bullet point expresses the fact that unitary evolution typically leads to thermal 
behavior, with an infinite temperature density matrix for this specific case. The second 
point expresses the deviations from thermality, showing that the error size when A ~ 0(1) 
is exponentially suppressed. This result has been used repeatedly in the past to argue 
against the Hawking information paradox [2], since the available computations of correla¬ 
tion functions during black hole evaporation are insensitive to the expected corrections, of 
O(e ^ BH ) if we trust the previous random approximations. 

We remark here, see [mug for recent reviews and references therein, that the intuition 
coming from the random unitaries/states framework has had a strong impact in the context 
of black physics and quantum thermalization, given a plethora of new concepts, problems 
and solutions. But we believe that the framework deserves further development, since its 
present status has also deficiencies. Firstly, the random unitary framework has not been 
formulated generically in the context of QFT. Secondly, the random unitary evolution 
defined in [El does not conserve any charge, such as energy, momentum, electric charge, 
etc. Besides, the framework just provides a Gibbs distribution with infinite temperature. 
It is then interesting to ask the following questions: 


• Can the framework be generalized systematically so as to include conservation of any 
given charge? 

• Can it be formalized in a QFT context? 

5 See jl2| . and references therein, for an analogous but somewhat different approach coming from quan¬ 
tum information theory. 
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• Can non-infinite effective temperatures emerge in the framework? 

• Do the deviations from thermality differ from the previous case? 

In this article we develop a framework to solve these questions, and which has interesting 
different potential applications. We summarize the results herej^j 

1.2 Summary of results 

The natural answer for the first question is not to choose the Haar ensemble of random 
unitaries/states, but to choose the appropriate ensemble of unitary matrices conserving any 
given set of charges. The problem with this approach is to construct such an ensemble, 
and we do not see a systematic way of doing so. In this article we will use the Haar ensem¬ 
ble in an appropriate subspace of the Hilbert space, defined by the given initial charges. 
This procedure can also be seen as implictly defining an esemble of random unitaries dif¬ 
ferent from the Haar ensemble in the full Hilbert space. Within this framework we will 
provide generic formulae for the typical reduced density matrix and entanglement entropy 
of any given subsystem and for any given set of conserved charges. These objects turn 
out to be written generically in terms of ‘constrained multiplicities’, number of different 
configurations with a given set of constraints. The results connect in a direct and transpar¬ 
ent way modern approaches of quantum information theory, concerning for example the 
computation of entanglement entropies, with more traditional frameworks of microstate 
counting. 

For the second question we argue that the natural way to generalize the framework to 
a QFT context is to use the Fock space formalism. Fock space provides an ideal setup 
to define the appropriate subspaces corresponding to the given set of conserved charges. 
Besides, it is also an ideal setup to define partitions/subsystems by tracing out different 
set of particles. One possibility is to integrate out all momentum modes but a given one, 
as done in [18] for the case of vacuum states. We can also integrate out positively charged 
particles, left movers of a CFT, etc. Interestingly, the formulas we find in the case of 
random states are generic for all types of subsystems. 

For the third and fourth questions we need to sacrifice generality of the analysis and 
choose concrete models in which to apply the developed formulas. We choose the simplest 
but non-trivial cases, in which many things can be computed analytically: a massless 
scalar field theory in 1+1 dimensions and its generalization to the case of N scalar fields. 
For a given initial total energy, we derive the typical density matrix and prove that in 
the thermodynamic limit it is equal to the thermal density matrix. It turns out that the 

6 Similar approaches to the one developed in this article have been explored in EH EH US [13 US]. 
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effective temperature is finite, and depends on the initial energy in the way expected from 
usual thermal considerations. In this way we are able to derive the equation of state of 
the QFT directly from random dynamics. The equation of state is seen as the typical 
macroscopic configuration of an underlying unitary quantum model. 

We also compute the deviations from thermality, obtaining an important insight. As 
commented before, for the so-called Page case, developed in m the typical density 
matrix is exactly the thermal distribution withouth the need of taking the thermodynamic 
limit, and the deviations of the typical entropy from the entropy of the typical reduced 
density matrix are exponentially suppressed. This means that deviations from thermality 
are exponentially supressed, see Appendix [A] for a complete discussion of all cases. On the 
other hand, as we will show, in the QFT context the typical density matrix and thererefore 
the entanglement entropy of the typical density matrix are not directly equal to their 
thermal counterparts, the difference being of 0(1/S). Therefore, to compute the leading 
deviations from thermality we just need to compute the entanglement entropy of the typical 
density matrix, a feature which will simplify matters in future developments. From other 
perspective, our results pinpoint the microcanonical ensemble as the right ensemble to use 
when considering pure state unitary evolution. 

At any rate, the important upshot is that the errors expected in a scenario in which some 
charges are conserved are much bigger than the ones commented before, being of 0(1/S) 
instead of 0(e~ s ). The applications of these results to conceptual problems involving black 
holes is still unclear and deserves further development. At any rate it is interesting that: 

• The corrections affect directly the diagonal entries of the density matrix, and therefore 
should be seen in observables such as mode occupation numbers. They are not non- 
perturbative, and should be visible in the first 0(1/S) corrections. 

• For any QFT, deviations from thermality can be computed exactly within our frame¬ 
work, and it would be interesting if this can guide the computations of these correc¬ 
tions directly from black hole geometric considerations. 

At the end of the article we speculate on potential applications of this framework to 
the relation between geometry and entanglement da ed]. In particular, our framework 
can compute entanglement entropies for any desired subsystem, and it is imaginable that 
we could get important insights into the emergence of near horizon regions through en¬ 
tanglement in the held theory. In previous literature there has been problems with the 
‘integration out’ procedure of high momentum modes, or the so-called ‘holographic renor¬ 
malization’, when trying to get to the near horizon region of black hole backgrounds, and 
therefore only the asymptotic region is well understood. It seems that our framework may 
overcome such difficulties in the specific context of random evolution, and might teach 
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us important lessons about event horizons, stretched horizons and near horizon regions of 
black holes. 


2 Reduced density matrices in Fock space and multi¬ 
plicities counting 

Our initial quantum state will be a pure quantum state within a Fock space, which is the 
natural basis of QFT and certain spin systems: 


|^aj ^bi } I®) 


( 2 , 1 ) 


The notation works as follows. Labels i a ,ib , ■ ■ ■ denote the number of particles in a given 
representation. Labels a, b, ■ ■ ■ specify the representation. For example, in a certain theory 
we could have a = p,s, q r , where p is the momentum, s the spin and q r some charges 


labelled by r. As usual we assume that the Fock states in (2.1) are attached with a total 


conserved energy Et and total conserved charges Q r T of the free theory, where r label 
different possible charges. Finally, when there is no need to write all the previous labels we 
will use greek letters a, f3 to simplify notation. They are just natural numbers running over 
the state^J For example, normalization of the basis vector is expressed as (a\/3) = 5 aj p. 

We now follow the logic described in the introduction. Consider that the full interacting 
Hamiltonian of the system H commutes with the total charges, so that E -r and Q ^ are 
conserved. If we begin with an initial Fock space state |T in ) = | (3), with definite values of 
Et and Q the various charge conservation laws force the unitary evolution of the initial 
state to be of the following form: 


<a=Q(F7t,Q'j’) 

!*(<)) = u(t)\p) = Y, • 


( 2 , 2 ) 


a =1 


where the sum does not run over all Fock states (2.1), but only for those |ck) with definite 
values of E t and Q r T . By definition there are Qt ) = of those, i.e the 

microcanonical number of states for a given total energy and total charges. 

As explained in the previous section, a useful approximation when considering chaotic 
and highly entropic sectors of the theory is to consider U{t > f relax ) — ^ ran dom- In this 
article, with the objective of extending the random unitary framework to systems with 
conserved charges, we will use Haar randomness in the previous subspace, which is indeed 
defined by the given set of charges. Up to second order corrections (doubly exponentially 


7 We suppose we have a discrete theory, probably with an infinite, but countable number of states such 
as in QFT. The examples we work out in the article fall in this class, since they are QFT defined on 
bounded domains. 





suppressed in the entropy of the given sector) this approximation can be defined by the 
following averaging procedure^ 


M = 0 = A5 afj ■ 


Imposing average normalization of the state we obtain 


(2.3) 


[(^out|^out>]=Afi(£ T ,QT) = l- 


(2.4) 


We want to remark that the ensemble of random states will always be the same, the previous 
gaussian ensemble of random states, which can be obtained by an analogue ensemble of 
random unitaries with random gaussian matrix entries. What will change is the sector 
of states in which gaussian randomness is considered. The sector of states will be the 
mathematical object which needs to be analyzed in several ways. From the perspective 
of the full Fock space, the relations (2.3) implicitly define the action of an ensemble of 
random unitaries, preserving any given number of conserved charges, on some state of the 
Hilbert space. 

With (2.3) we can now ask questions about typical properties of the held theoretic state 
| T out ). For example, the average of the global density matrix is given by 


oi=Q(E r Y ,Q^) 

[pout] = [|^out)(^out|] = ^ nr , (2.5) 

which is just the microcanonical density matrix for a theory with a generic number of con¬ 
served charges, obtained here in a straightforward manner by applying random dynamics 
in the appropriate subspace. Random dynamics seem to pinpoint the maximally mixed 
microcanonical state, instead of the canonical ensemble. 

The next natural step, given the structure provided by Fock space, is to choose one 
particle type, say a — p,s, (f, and integrate out all the others. Notice that the case of mo¬ 
mentum space entanglement, analyzed in jTS] for vacuum states, is going to be seen here as 
a special case, arising when we consider all types of particles within a specified momentum 
shell. But otherwise the Fock space formalism provides naturally more possibilities, such 
as the entanglement between positively charged particles and the rest, particles with a 
given spin and the rest, left and right movers of a CFT, etc. All possible types of sets and 
reduced dynamics within the Fock basis of the Hilbert space are constructed by joining the 
appropriate set of a, b, ■ ■ ■ conforming the type of subsystem we want to study. We thus 
begin with the simplest one, given by one particle type a. We first write the state in the 

8 In Appendix |a] we show how to derive Page’s formula for the average entanglement entropy 5] with 
these simple Gaussian relations. 
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following form 


( 2 . 6 ) 


|*out> ^ ^ ^iq,a\^ai Q^) i 

ia 

so that now a represents all other particle types different from a, but the sum still runs 
over the previously specified subspace. The reduced density matrix is 


Pa = 


^ ^ (o^| fl/put) (fl/put |Q() — ^ia,oi^i a) o)Va) {}a\ ■ 


(2.7) 


This is directly a diagonal density matrix, without the need of averaging. Indeed, for a 
subsystem with one particle type a, it just happens that for two states | i a , a) and \j a , a) 
with the same energy and charge, they must obey i a = j a ■ It will no longer be true for 
subsystems with more than one particle type. Using (2.3) we can obtain the typical reduced 
density matrix 


w = E 


V- (lUx-Cr.',) ■ 


H n(£?r,gy 






( 2 . 8 ) 


ia= 0 


where Q(E T , Q^, i a ) is the number of states with energies E T , charges Ql r , and a fixed 
number i a of particles of type a. The maximum value of the number of particles of type a 
is denoted by i™ ax and needs to be determined case by case. We discuss this in the example 
presented in the next section. 

With the previous density matrix, the average entanglement entropy of a particle of a 
given type a, in the fixed E t and Q r T sector is given by: 


[-S'e(Po)] 


;max 

L a 

p (^T, <?T, ia) log P{E T , Q r T , i a ) . 

ia =0 


(2.9) 


Being rigorous this is the entanglement entropy of the average density matrix, usually de¬ 
noted by Se ( [p a ]). The full average entropy is computed in Appendix [Aj where we show 
that the difference between the two is exponentially suppressed in the microcanonical en¬ 
tropy. The entropy of the typical density matrix will be enough for us, since we are not 
interested in the deviations between the typical entropy and the entropy of the typical 
density matrix, but instead, on the deviations from thermality. In the usual case of [5], 
typicality and thermality turn out to be exactly the same, and therefore deviations be¬ 
tween thermal entropy and average entanglement entropy are equal to deviations between 
average entanglement entropy and entanglement entropy of the average (see Appendix [A] 
for the complete discussion). We will show that when conserved charges are taken into 
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account this is not longer true. The leading deviations from thermal entropy are those 
already accounted for by the entropy of the average density matrix, a feature which sim¬ 
plifies present and further developments. At any rate, if needed, the full computations are 
described in Appendix |A| 

Relations ( |2.8 ) and (2.9) are generic formulas which turn out to be written just in terms 
of ‘constrained multiplicities’, the number of different states with a given set of constraints. 
Notice that the result for the entanglement entropy is finite, and no divergences occur, even 
taking into account that we are dealing with quantum field theories. This is obviously 
because we are carefully applying conservation of charges. Although the multiplicities 
might be difficult to compute in general, there might be a class of theories in which they 
can be computed, and explicit connections with thermal density matrices (in particular 
generalized Gibbs ensembles) might be established. We do this in the following section for 


the case of certain scalar held theories. At the same time the generality of (2.8) and (2.9) 


suggests there might be a generic way to prove that the probabilities P(E/r, QtGg) are 
well approximated by the Generalized Gibbs Ensemble. Also it would be interesting to 
compute this multiplicities for integrable theories. In those theories Generalized Gibbs 
Ensembles are expected to govern the dynamics but deviations from them might be bigger. 
This might aliviate several problems encountered in previous literature, see [7j. We leave 
these interesting paths for future work. 

Generalizing the procedure to include any desired subset of particles is straightforward. 
We first form the set of a, 6, • • • , c in which we are interested, and write (Tout) as: 


l^out) 'y ^ db a ,4,-" j c:CX \i a Afe; , i c ; o;) , 


( 2 . 10 ) 


^ct j*'' 


where a now labels all particles types different from a, b, ■ ■ ■ , c. The reduced density matrix 
then reads: 


Pa,b, -,c y ^ (o. | fl/put) (if^out | q) 


( 2 . 11 ) 


N (E J V' 2.. 

i,i' \ a / 


• \ / •/ •/ •/ i 

? ^c) I 5 


where a labels the set of particles which are traced out. As opposed to the previous one- 
particle case, this is not a diagonal density matrix^] However, taking the average we find 
the following generic formula: 

r i Qti ®6i ’ • ‘ Ac) |. • w • . . | 

[Pa,b,---,c\ / c , / rn \ PaGb) j k)(*a) ' Gc| 


9 Although it is not diagonal, it has a nice block diagonal structure, as used and described in Appendix [a] 
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( 2 . 12 ) 


^ ^ P(E T j Q'j'i iai ibi > ic) \ ^at ^bi Ac) (iai ^bi Ac | i 
i 

which is diagonal. The average probabilities P(E Tj Ql r , i' a , i' b , ■ ■ ■ ,i' c ) are ready to be com¬ 
pared with thermal expectations, i.e with probabilities coming from Gibbs distributions. 
We will discuss examples in the next section. 

The average entanglement entropy - i.e. the entanglement entropy of the average density 
matrix (see Appendix [AJ- is finally 


[SeG^a- ,c)] = - P (^T, Qt Aa, ib, ■ ■ ■ , A) log P(E t , Q r T , i a , i b , ■ ■ ■ , i c ) . (2.13) 


We concldue that the generic equations (2.12) and (2.13) apply to any theory with a Fock 


basis structure, such as a QFT or certain spin systems. In particular it is applicable to 
Holographic Field theories. It seems a rigorous framework to study entanglement between 
infrared and ultraviolet domains at finite temperature. An exciting possibility is that we 
could potentially extract physics from the near horizon regions of quantum black holes 
from the structure of entanglement in Fock space of the field theory. At any rate it is 
expected to give more insights into the connections between entanglement and quantum 
gravity PH no, an exciting direction which we leave for future work. 


3 Examples: massless scalar fields in two dimensions 

In this section we analyze in detail two specific examples. The first example is a massless 
scalar field theory in 1+1 dimensions defined on a finite line segment. The second one is its 
generalization to the case of N scalar fields. For both cases we will show the emergence of 
Gibbs distributions as effective states for reduced subsystems. The procedure also allows 
the computation of deviations from precise thermality. Finally, entanglement entropy for 
single modes as a function of the momentum, and entanglement between high and low 
energy momentum modes for a given energy cutoff will be discussed. We will end with the 
analogue of Page curve (5j for the case at hand. This enlarges the program spelled out in 
[18] to the case of random states. 

3.1 Massless real scalar field on a finite segment 

In this case the energy/momentum dispersion relation for excitations over the vacuum, 
together with the quantization condition reads 

7m 

E n = p n = — , (3.1) 
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for a segment of length L and n = 1, 2, • • ■. The Fock space is spanned by vectors of the 
type: 

I®) = |fy=l; fy=2j ' ) |fy) fy) ' ) j (3.2) 

where i n is the number of particles with momentum p n , and a is just a natural number 
running over all the infinite but countable eigenstates, used here to simplify notation. 
These are eigenstates of the free Hamiltonian with Dirichlet boundary conditions. If the 
true Hamiltonian conserves the total energy, which for a Fock state reads 

71 =Umax Umax 

E T = ^ inPn = -jr Y ^ , (3-3) 

n =1 n =1 

where n max is defined by p nmax = Et, then any initial state with definite total energy Ify 
will evolve towards states of the form 


V I'out) = d'afy) , 

a=l 


(3.4) 


where the |a) are all the states belonging to (3.2) with total energy Ffy. By definition these 
are Q(Et) = e s ^ ET \ where S(Et) is the microcanonical entropy at energy Et- Noticing 
that (|3.3|) can be written as 

Y i n n , (3.5) 


-E t = 
7 r 


n=l 


we conclude that the number of states with a given Et is equal to the number of different 
partitions p(n) of the natural number n max = -Et, which is given asymptotically for large 

«max by 




fl(E T ) - p(—Et) -> r- -< 

7T A\/3EtL 

The global entropy can now easily be extracted and for large n max we find 


(3.6) 


S = 



(3.7) 


So, in the limit Et —> oo, this agrees with the Cardy formula for a CFT with central charge 
c = 1 and L 0 — ^ = n max . If the interacting theory would be a conformal theory on a 
line segment of length L , the conformal dimension of the pure state T out ) would then be 
A Umax' 


Now we can directly apply the generic formulas derived in the previous section. The 
approximation U(t > t re i ax ) — Random is operationally defined by (2.3), which we repeat 
here to emphasize that it does not change from one case to another: 


M = 0 IK*?] = A <W ■ 


(3.8) 


13 












Imposing average normalization, we obtain [('bout|'bout)] = AQ(Et) = 1. The average of 
the global density matrix is given by 

a=fl(E T ) 

[Pout] = [|'bout) ('bout |] = ’ ( 3 ‘ 9 ) 

which is just the microcanonical density matrix, obtained here in a straightforward manner 
by applying random dynamics in the appropriate subspace. 


3.1.1 Entanglement of a single momentum cell 


We now study the entanglement entropy in a pure state of a single particle specie. In our 
example, choosing one particle type just amounts to choosing a definite momentum p n for 
a fixed n. The momentum cell can still be multiply occupied, so the subsystem has a finite 
dimension determined by the size of the cell, i.e. i™ ax + 1, where the maximum occupation 
number is given by i™ ax = [E T /p n \, the integer closest to E T /p n from below. Writing the 
state in the form 

l^out) = , (3.10) 

in tOL 

where a now represents all other particles types different from p n , we arrive at 

.•max 

L n 

Pn = £(«l* OU t)(*out|a) = ^ ] 'bj W]Q 'b iin , Q ,)|i n )(i w | . (3.11) 

a in=0 a 


This is directly a diagonal density matrix, due to energy conservation. The average density 
matrix is 


w = E 


1=0 


p, in) |. v /. 


•max 

L n 

^ ^ ^ in)\in) {^n\ •> 

in— 0 


(3.12) 


where Vt{E^,i n ) is the number of states with total energy E t and i n particles with mo¬ 
mentum p n . Therefore, the average entanglement entropy of momentum cell n is given by 


[Se(P„) 1 = - E P ( P T, >») lQ g P ( p T. in) ■ 


(3.13) 


in —0 


Formulas (3.12) and (3.13) are the analogues of (2.8) and(2.9) for the case at hand. To 
compute these quantitities we need to find f2(i?T,in)- In this example this number can be 
analytically computed. The procedure is explained in detail in Appendix [Bj Here we just 
quote the result: 


H(E t , i n ) = p ( ^ E t - i n n j - p (-^E t - n(i n + 1) 


(3.14) 
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where p{n ) is the number of partitions of the number n. 

Now we are ready to make the first crosscheck of our computations, which constitutes 
one of the main results of the article. The typical probability of finding i n particles with 
momentum p n is given by 


P(E T ,i n ) 


, in) 

n(E T ) 


p(^E T - i n n ) - - n(i n + 1 )) 

P(^ E t ) 


(3.15) 


It is easy to check that all probabilities add up to unity. The leading term when i n p n £T, 
or equivalently i n <C i™ ax , is given by 


P(3 


Tj *n) 6 


~Pn (*n + l) 


(3.16) 


This is a thermodynamic limit. Indeed, i n p n Ft means that the energy in momentum 
cell n is small compared to the total energy, or equivalently, momentum cell n is not well 
occupied. This implies that most of the energy is sitting in the other momentum cells, and 
that these cells form a heat bath for momentum cell n. Another way of saying this, is that 
this limit is a good approximation for low values of n. For higher values of n, the energy 
in this momentum cell will typically be too large, or = \_n max /n\ will be too small for 
the thermodynamic limit to be a good approximation. If we would have used the usual 
Gibbs ensemble this probability would be given by 

p flinPn 

P(T, i n ) = ™-- 5 — = e- e -' 5 »-‘ i " +1 > . (3.17) 

E<m =0 6 

So in the limit i n p n <C E T , of for CFT’s, large conformal dimensions A = n max , the pure 
state is typically seen by the momentum cell as a thermal bath at a temperature given by 


T = 



(3.18) 


We thus see that it is possible to derive generic Gibbs ensembles, with any desired effective 
temperature, using random dynamics in the approppriate subspace. Notice that (3.18) 
implies S(Et)T(Et) = 2E t on the one hand, and on the other hand there is the general 
relation S(Et)T(Et) = E t + VV. Combining the two, we can determine the pressure 
density: 

V = E t /V = e . (3.19) 


This is the same equation of state as for a two-dimensional CFT. In this way we expect 
to recover the known relation e = (d — 1)V, valid for d-dimensional CFT’s, using just 
random dynamics. More generically it should be possible to derive any equation of state 
by choosing the appropriate field theory and using the same procedure. The equation of 
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state of a fluid system, at least for this case, is the typical macroscopic configuration of 
the true evolving pure quantum state. 

The next step is to compute the deviations form thermality. The next to leading term 
in the previous i n p n <C E T expansion is given by 


P(E T ,i n ) = P(T,l n ) + + l) e -frn(in + 1) = 

j Ex E t 


= P(T, i n ) + error . (3.20) 


To proof this, we only needed the Hardy-Ramanujan asymptotic formula for the number 
of partitions, see (B.12) in Appendix |Bj 

The error we find here is thus much bigger than the one is usually expected. Within 
the usual random unitary formalism, the errors are typically of 0(e~ s ), so exponentially 
suppressed in the entropy. In our case, if we take a momentum cell with energy i n p n ~ T, we 
are finding errors of 0(T/E) ~ 0(1/S). This seems just to be due to energy conservation, 
which is explicitly ensured in our formalism. That we find such large errors may have 
implications for black hole physics and bulk locality in AdS/CFT, see [TU 12] and references 
therein, and we address this issue further in the next subsection when we discuss the large 
N limit. 


Finally, the average entanglement entropy of the momentum mode is given b} 10 

.•max 

L n 

[S E (n, ivj] = - ^2 P( e t, •») log P(E t , i„) , 


(3.21) 


—0 


and we remind that i™ ax = |_ E T /p n \ = [n max /n\. The dominant contributions to this sum 
come from the low occupation numbers, where the probabilities are approximately thermal. 
For high occupation numbers, i n ~ i™ ax , the probabilities are exponentially suppressed in 


the entropy, as one can easily compute from (3.15) 


At present we have not found a way to evaluate the sum algebraically in closed form 
as a function of n, though one can evaluate the sums explicitly for any given n and n max , 
e.g. on Mathematica. However, for high values of n this is possible. For example, in the 
case of n — n max , the highest momentum possible, there are only two terms in the sum in 
the entanglement entropy, i n — 0 and i n — 1 = i“ iax , and the result is 


[S E (n n 


J] = - 1- 




log 1 - 


log 


p( n max )) p{ ^max ) p(n 

max) 


Se 


-s 


(3.22) 


where in the second line, we took the leading term in the limit n = n max —)■ oo. The 


relation (3.22) shows that there are extremely small entanglement entropies in random QFT 


10 In the sum it might happen that some probabilities are zero. These terms simply do not contribute 
to the sum, since lim x _>o x log x = 0. 
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states 11 For conformal field theories, and in the context of the AdS/CFT correspondence, 
we do not expect these entanglement entropies to be captured by some geometric quantities 
in the bulk. The proposal for deriving entanglement entropy of CFT’s holographically, 
developed in [H5], is expected to capture entanglement entropies with a minimum size of 
(9(1), since this would correspond to surfaces of Planckian size. Entanglement entropies 
of 0(Se~ s ) are clearly of non-perturbative nature from the point of view of AdS/CFT 
dualities m 

On the other hand, the thermal entropy of a single mode is based on the Gibbs prob¬ 


ability distribution (3.17). The result is well-known and can be directly computed using 
P(E' r, i n ) and Shannon’s expression for the entropy. It reads 


3-PPn 


Sp{n) = pp n j ^ e -p Pn ~ lo §( 1 - e /3P ") ’ 


(3.23) 


with inverse temperature /3 = l/T given by (3.18), and we remind that p n = irn/L , such 
that (3p n = nn /v/6n max , which is independent of the size L. The function (3.23) is a 
monotonically decaying function. 

The two functions, the entanglement entropy and the thermal entropy, are plotted 
against the momentum n in Figure [lj The dependence against n provides some short of 
entanglement ‘running’ on the energy scale of a thermal-like state. 



Figure 1: Average entanglement entropy S-e(ii) (left) of a single momentum cell n and thermal entropy 
Sp(n) (right) seen by the momentum cell, for n max = 1500 and n varying between 1 and 200 (horizontal 
axis). For this large value of n max , the agreement is very good. 


It would be interesting to End an asymptotic formula for the entanglement entropy for 
large and generic momentum p n . This would make easier the comparison with the ther¬ 
mal approximation. Numerical analysis, as can be seen from Figure [2j shows that for large 
values of n max , the entanglement entropy approaches the thermal entropy, but the correc¬ 
tions to the thermal result are not exponentially suppressed in the entropy, and go like 

11 Notice that this is not the case when considering the common Page case J5] . Indeed, in that case the 
minimum entanglement entropies are of 0(1) in the thermodynamic limit. 
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1 /ydz max ~ 1/A. This was obviously expected from the corresponding errors in the prob¬ 
abilities themselves. At any rate, notice that since corrections die in the thermodynamic 
limit, the expectations coming from the simple analytical thermal expression (3.23) might 
be ready to compare with geometric dual formulations of the QFT, a very interestring 
direction to explore. 



Figure 2: Difference of entanglement and thermal entropy, now as a function of n max , and for fixed n here 
taken to be n = 2. A best fit shows that the fall-off of this graph goes approximately like l/y / h ma x ss 1/S, 
and not like e~ s . As explained earlier in the text, this only holds for low values of n. For large values of 
n, the differences become larger. 


3.1.2 Entanglement entropy of low-energy degrees of freedom and Page curve 


As described before we can also consider any subset of particles we wish. In this case a natu¬ 
ral property to analyze is the entanglement between high-momentum and low-momentum 
modes for a given threshold p n . More explicitly we can integrate out all particles with 
momentum higher than p n . The reduced density matrix is then given by 


P 1,2, 






Vi/* 

V ^ , 


pi, * 2 , 


A(A, ^2, 


i' I 

? L n I 


(3.24) 


i,i' ol 

where a labels the set of particles with momentum greater than p n . As in the previous 
section, this is not a diagonal density matrix. Taking the average we find 


r i _ n(^,*l,*2, • • • An) | . . . w. . . | 

[pl, 2, — / y Q(E ) rlj *2) ■ ■ ■ ) l n)( l l) l 2,’" An | — 

^ ] ^~ > (A'tAiA 2, An) |*lA2,‘'' An)(hA 2, ' ' ' An| • 


(3.25) 


The average probabilities P(E T ,i 1: i 2 , ■ ■ ■ ,i n ) can be computed using the same previous 
method of constrained partitions. This is explained in detail in Appendix |B} The result is 
given by (B.17) divided by the total number of partitions p(-^Et): 


P(EtAiA 2 ,-- - An) = 


P 


Q(£ t ) 2 tt 


-E t - ^ ikk ) + 


k=1 


18 










+ 


+ 


i>( 

;=i \ 

D>( 

l<r \ 

(-1 Tp 


L 

-Eh-J2i,k- 

k=1 

L ^ , 

k=l 

\ fc=i 


+ 



+ ••• 



(3.26) 


This is a quite complicated expression for generic though one can evaluate 

these sums using Mathematica rather easily. The situation becomes more cumbersome 
when evaluating the entanglement entropy, since we now have to perform additional sums 
over multiple occupation numbers. In practice this turns out to be rather hard using 
Mathematica, and we only succeeded to evaluate the entanglement entropy of a system 
of up to the first eight modes. However, a simplification occurs whenever Y^k^kPk •C E- r, 
namely when the subsystem has much smaller energy then the total energy. Notice that if 
r <C n, then we can approximate the partitions 


p(n — r) ~ p{n)e ^ , 


(3.27) 


which we can use in the probabilities (|3.26). Indeed, if )Cfc=i ' l kPk "C E r holds, we can 


apply (3.27) to all terms in (3.26), so that we obtain 

g — P Ek= 1 ikPk 


(3.28) 


— g ^^2k=l ikPk 


]^[(1 ~ e~ 0Pl ) = P(T,ii, i 2 , • • • , i n ) , 


i=i 


where P(T,ii,i 2 , ■ ■ ■ ,i n ) is the probability given by the Gibbs ensemble. One can again 


compute the subleading corrections, using (B.ll) for each term in (3.28) and find power 


law suppressed terms in the energy of the subsystem divided by the total energy. So the 
error analysis is similar to the case of a single momentum cell, and the leading corrections 
scale like 1/S. Hence for any subsystem, up to computable corrections, the entanglement 
entropy associated to pig,... is just the sum of thermal entropies of each momentum mode. 
In this way we can derive the analogous curve in a QFT setup to the so-called Page’s curve 
[SJ. Obviously, if we express the entanglement entropy as a function of the corresponding 
thermal entropy we obtain directly Page’s curve, with somewhat different deviations from 
thermality. But in this context it is more interesting to paint the entanglement entropy as 
a function of the energy scale used to divide the high-energy modes from the low-energy 
ones, given by p n . Hence we define two complementary subsystems A and B in momentum 
space, as depicted in Figure [3j 
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A B 


1 


^ ^max 


Figure 3: Two complementary subsystems A and B in Fock space. The entanglement entropy of these 
subsystems is computed below. 


Up to subleading corrections in the thermodynamic limit, we can use the thermal 
entropy for the low-energy modes, which becomes simply the sum of the thermal entropy 
of the individual modes 

SaH = [P pk 1 1 e -p Pk ~ M 1 ~ e ~^ Pk ) ) • ( 3 - 29 ) 

k =1 k =1 ' ' 

The plot of the thermal entropy S A is given in Figure [4} It is hard to perform the sum in 



Figure 4: Thermal entropy (vertical axis) of a subsystem A containing all momentum modes up to a 
given n, as a function of n (horizontal axis). In the figure we took 1 < n < 200 and n max = 1500. The 
curve asymptotes, for n —» n max , to 96.22, which is reasonably close to the total entropy of the system 


given by (3.7), which for n max = 1500 yields S ~ 99.35. The discrepancy gets smaller for larger n n 


The entanglement entropy approximates the thermal entropy well whenever the energy in subsystem A is 
small compared to the total energy, so for small values of n. Deviations then go like 1/S. 


(3.29) analytically, except in the continuum limit where the modes /3p n become infinites¬ 
imally spaced. Since (3Ap = tt/\/Q n ma ^, the continuum limit is n max -A oo, or L —* oo 
keeping temperature fixed. In this limit we can approximate the sum by the integral 

i rPp™ 


S A (n) = 


/3Ap 


da; x 


S 

2 


'o 

7m 2 


1 


65 (3Ap 


x ~ log(l - e~ x ) j + 0(1/L) 

(Li 2 ( 1 - e^ Pn ) + Li 2 (e- ppn )) + 0(1/L) , 


(3.30) 
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where we used the indefinite integral 


dx 


x- 


log(l 


x 


= -y-r*2(l 


e x ) - Lh(e~*) , 


(3.31) 


and S is the total entropy (3.7). As a crosscheck one can verify that in the limit n —> oo, 
one obtains the total entropy. This can either be seen from the asymptotic expansion of 
the dilogarithm, or by directly doing the integral 


/ 7°° / e~ x \ 

S = ~pl ix ~ log(1 " e ~’ ) ) ■ (3 ' 32) 

We can also determine the thermal entropy of the complementary system, which we denote 
by 5, namely the system of momentum modes that lie between a given n and n max . This 
is given by 

77 max 71 m ax / _ fipfc \ 

S B (n) = E ft#) = E (ft* t Z ~ Ml - c-‘ ,n ) ) • (3-33) 

k=n k=n ^ ' 

Applying the previous continuum limit we can obtain also analytical formulas for the 
entropy of B. The plot of the thermal entropy of subsystem B is given in Figure (J5]). The 



Figure 5: Thermal entropy Ssin) (vertical axis) of the complementary subsystem B containing all 
momentum modes between n and n max , as a function of n (horizontal axis). In the figure we took 
1 < n < 200 and n max = 1500, just as in Figure [4] The curve now asymptotes for n —>• 1, again to 96.22, 
which is reasonably close to the total entropy of the system given by (3.7), which for n max = 1500 yields 
S — 99.35. The discrepancy gets smaller for larger n max . The thermal entropy of system B is reduced 
to half at mode number around n ~ 21 again. For larger values n 21, the thermal entropy is a good 
approximation for the entanglement entropy, with corrections of order 1/S. 


thermal entropies associated to A and B have different functional structures as we vary the 
subsystem sizes. Given that the entanglement entropy of A should be equal to that of B 
it might naively seem there is a contradiction here. But indeed there is no contradiction, 
since the thermal entropy is a good approximation for the entanglement entropy of A for 
values up to n — 21 (half the entropy), and for the entanglement entropy of B for larger 
values of n. 
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Therefore, what is meaningful here is the critical momentum mode p n for which the 
thermal entropy of the reduced subsystem is precisely half of the total entropy. This 


happens when the corrections to the first term in (3.30) cancel between each other. Re¬ 


markably, in the large L limit, there is an exact solution of this, due to the identity of 
dilogarithms 

Li 2 ( 1 — x) + Li 2 { 1 - aT 1 ) = -^(logx) 2 , (3.34) 


which we will use for x = 2 and /3p n = log 2 in (3.30). For these choices, all other terms 


cancel except for the first one. Hence the critical momentum dividing the QFT in two halfs 
maximally entangled with each other is given by 


p crit = T log 2, 


(3.35) 


giving a different interpretation of temperature in the QFT. From this perspective, the tem¬ 
perature T provides the energy scale wich divides the QFT into two equal parts maximally 
entangled with each other, with a entanglement entropy equal to S/2. 

In terms of the mode numbers, (3.35) can be written as 


n cnt = V6log2 —— _ o.54^n max . (3.36) 

IT 

This can be checked by looking at Figure a for which we have n crit = 0.54 a/ 1500 ~ 21. 
For this value of n, we have SU(21) = 47.33, whereas S/2 = 49.67, so indeed around half 
the total entropy. 

The Page curve, as a function of the energy scale dividing high and low energy degrees 
of feedom, can now easily be found by combining the two curves in Figures [4] and [5] They 
intersect exactly at n cnt ~ 0.54^/f hTT , as depicted in Figure [6j 



Figure 6: The ‘Page curve’, i.e the average entanglement entropy of reduced subsystems in random 
states, as a function of the energy scale dividing high and low energy degrees of freedom. The turning 
point, where the entanglement entropy takes its maximum value S/2, occurs for an energy scale of 0(T ), 
providing a different interpretation of temperature in the field theory. 
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3.2 The large N limit 


As commented in the introduction, part of the reasons to generalize the random unitary 
framework to the context of quantum held theory is to get closer to the physics of black 
holc^J In particular our formalism can be useful to describe black holes in anti de Sitter 
spacetimes via AdS/CFT [51J- In these theories we have a large gauge group, with a 
corresponding large number of held species. In AdS/CFT this number of held species is 
counted by the central charge c of the CFT, which is taken to be large to have a smooth 
gravitational dual background. With this in mind we would like to repeat the previous 
exercise for the case of N scalar helds, in the large N limit. 

For each scalar held, labelled by a = 1, • • • , N, we have the same energy and momentum 
dispersion relation for excitations over the vacuum 


E n a — Pn a — 


Tm a 


(3.37) 


for a segment of length L and n a = 1,2, • • •. The Fock space is spanned by states of the 
type 

|®) = \^n a i ^n b t ' ) ) (3.38) 

where i Ua is the number of particles of the scalar held a with momentum p Ua , and a is again 
a natural number running over the infinite but countable set of eigenstates, used here to 
simplify notation. Conservation of total energy 


CL —iV Tla — Umax 

e t = in - Vn - 

a— 1 n a =1 


a=N n=n max 

£ 


a =1 n a =1 


,n n 


(3.39) 


where n max is as in the one-held case p nmax = At (so n max = LEnc/n), implies that any 
initial state with def ini te total energy E t will evolve towards states of the form 


QL=£l(E r Y , AT) 

l^out) = V a \a) , 

a=l 


(3.40) 


where |a) run over all the states belonging to (3.38) with total energy Et- By definition 
these are Q(E^,N) = e s<yET,N \ where S(E^,N) is the microcanonical entropy at energy 
Et of the theory with N scalar helds. In Appendix |B| this number is computed, and in the 
limit Et Nn/L or equivalently, n max 3> N, we obtain 


N 


Q(E t ,JV)=(p(^)] -> 


N 


4'Umax \/3 


N 




(3.41) 


12 This formalism might be also interesting in the context of integrable theories and Generalized Gibbs 
Ensembles, see j7|. 
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consistent with the Cardy formula for the microcanonical entropy: 


S(Er,N) ~ J^LEtN = ZnJ^cA , 


A = L 0 - — = n max c — N . (3.42) 


Now we can repeat the same generic procedure described in Section [2] In this case we will 
apply the generic formulas for the average density matrix and entanglement entropy (2.12) 
and (2.13) directly. The average of the global density matrix is given by 


[pout] — [|^ r out)(^ r out|] — h2(_E T 7V)l a )^l 


which is again the microcanonical density matrix. 

The typical density matrix associated to a de fini te momentum p na is 


(3.43) 


[ftj = E 


—0 


fl(E T , A, L . 
n(Er,N) 


in a )(in J= N,i na )\i na ){i na \ , (3.44) 


inn —0 


where N,i na ) is the number of states with total energy Et and i Ua particles with 

momentum p na , and i™ ax = \_E^/p na \. Therefore, the average entanglement entropy of a 
particle with momentum p Ua is given by 






P(E t , N, in a ) log P(E t , N, i r 


,=o 


(3.45) 


Formulas (3.44) and (3.45) are again the analogues of (2.8) and (2.9) for the case at 
hand. To analyze the previous formulas we need to find the constrained multiplicities 
0 (£'t, N, i na )- This is computed in Appendix |B| in the same limit E t Nn/L as before. 
The result reads 


r)( T? at ■ \ ( / max * n a n a\\ ( Amax n a(in a "b 1)\\ / Q 

n{E T ,N,i na ) = I p( -—- ) 1 - I p{ -—-) I , (3.46) 

where we remind that n max = LE^/n. Expanding the probabilities P(E^, N,i Ua ) for 
i na n a n max , the leading term is given by 


P(E T ,N,i na ) -> e~^ inaPna - e" 
which is just the Gibbs ensemble at temperature 

rp _ \ _ f 6-^T 

/ 3 V nNL 


(3.47) 


(3.48) 
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Notice that again TS = 2Et, so that the pressure density reads 


V = E t /V = e , 


(3.49) 


as expected for a two-dimensional conformal held theory. 

We can now compute the deviations form thermality, to check if there is some extra 
dependence on the central charge c = N of the theory. The next to leading term in the 
previous expansion is given by: 

P(E t , N, i n J = P(T, N, i n J + N l ^E e -^ Pn - - N Pna ^ na + ^ e ~^ ( ^ +1) . (3.50) 

Et Et 

The error we obtain is again much bigger than usually considered. For a typical momentum 
mode with energy of 0{T), the error^ T (T) is: 

NT N 1 


error^T) 


rv.; - rs./ 


E t S(E t , N ) TL ' 

In conformal held theory language, this translates into 


(3.51) 


error^(T) 


1 

TL 


(3.52) 


so for large conformal dimensions A compared to the central charge c, the errors are small. 
Equivalently, for high temperatures, the errors are small. 

There is a simple intuition behind this result, given the approximation in which we 
are doing the computations. As explained in Appendix [Bj to compute the microcanonical 
degeneracy LI(Et, N ) for the case of N helds, we need to perform a difficult sum over 
different partitions. To compute it, we use a saddle point approximation. The physics 
behind the saddle point approximation, expected to be valid for E t Ntt/L , is that 
the problem at hand is equivalent to N single scalar held theories, each one with a total 
energy E t = Et/N . In this approximation we can then use the results of the single held 


model with total energy given by E t- Formula (3.20) then provides an error of order 


iPrig 


Et 


= N ln ^ a , and indeed, in the limit E t Ntt/L, the error is small. 


We want to remark here that one should not believe (3.51) at all energies, since the 


previous formulas for the multiplicities are only valid in the limit Et Ntt/L, or equiva¬ 
lently, for high temperatures TL 1. Naively it might seem that thermal physics should 
not be valid for E t < Ntt/L, but this is not the case. For smaller energies it is difficult to 
do the exact computations. But to see that thermal dynamics is still a good approximation, 
and compute the deviations from it, we can consider the specihc case in which the total 
energy is given by Et = tt/L. The only Fock states with such an energy are the ones with 
one particle excited in the lowest momentum mode. There are N of such states, one per 
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field specie. Therefore, using the previous random machinery, the typical reduced density 
matrix of the lowest momentum mode of a single field is given by 


P(Pi) 


N- 1 
N 


|0pi)(0 P i| + 


1 

N 


lpi)(lpi 


(3.53) 


This is just a mixed state, with probability equal to for specie 1 to be in the vacuum 
and probability equal to A of having its lowest momentum mode occupied. Matching 
with the associated Gibbs ensemble is not as transparent as before. One possibility is to 
match the first probability exactly, from which we obtain (3 = ^ log N. With this effective 
temperature the other probability is given by A(1 — 1/A). We see that the deviation from 
thermality is of 0(1/N) in this specific case, so we conclude that thermal dynamics is still 
a good approximation to the expected unitary microscopic result in the large N limit. At 
any rate, it seems there might be a qualitative difference in the analysis when going from 
energies Ntt/L to E T < Ntt/L, or equivalently from high temperatures to below a 

critical temperature T crit ~ 1 /L. Assuming that our results would also be valid for a CFT 
on a circle, and the field theory at large N is dual to a gravitational bulk with AdS 3 radius 
proportional to L , then this would signal the Hawking-Page transition 1221, seen in the field 
theory in the size of deviations from thermality. This might have interesting implications 
in the context of AdS/CFT [23]. The regime we are working in corresponds to the black 
hole phase. Indeed, the black hole corresponds to a highly excited state in the CFT, so 
we require a large conformal dimension A. Furthermore, large central charge is required 
for the bulk theory to be weakly coupled, and finally, we require A > c or T > T crit to be 
well above the Hawking-Page transition point. 

Moving forward, the average entanglement entropy of the single momentum mode n a 
is given by 

-•max 

n a 

<SE(n a ,n max , N) = - P(E t , N,i Ua ) log P(E t , N,i n J . (3.54) 

As stated before, this is a complicated sum for which we do not have a definite analytic 
expression. For the case of n a = n max one can explicitly do the sum, since there are only 
two terms, i na = 0 and i Ua — 1. The result is 


Ee,N (^inax) = —II — 


p(/7 m ax • N) 

log 


l0g 0 P(™max, N) 


p(n 

max 5 N ) p(^max) AT) 


S( E, N) e - s ^ N) , (3.55) 


where in the second line we took the leading term in the limit n = n max —> oo. We see 
that including more field species does not spoil the appearance of exponentially suppressed 
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entanglement entropies in momentum space. We remark again that we find it very improb¬ 
able that geometric quantities, such as geodesic lengths, can capture such physics, since 
the minimum proper lengths are of 0(1) in Planck units. 

Other aspects of the entropy behavior do not change drastically either from the single 
field case. This is easily observed by noticing that the thermal entropy of a given momentum 
mode, whatever the field specie, is left basically unchanged by the inclusion of more field 
species. The only change comes again by the total energy associated to a given field, and 
therefore the effective temeprature associated to each mode. For N fields, the total energy 
given to one field is E T = E T /N, and therefore the temperature T = = J fyy 

decreases with N as well. The deviations are given again by = N ln °£" a ~ N/S, as 

the previous probabilities themselves. 

Integrating out all momentum modes higher than a certain critical momentum p n is 
also possible. The associated constrained multiplicity can be computed in the same way 


as the previous ones, and an analogous formula to (3.26) can be obtained. Therefore, one 


can approximate the typical reduced density matrix with the Gibbs ensemble in the same 
previous limit. We will not repeat the procedure here, since it is straightforward and those 
results are not affected by the inclusion of more fields. 

Besides, the analogous of Page curve in this QFT setup is not going to change either. 
To observe this, notice that up to subleading corrections in the thermodynamic limit, we 
can use the thermal entropy for the low-energy modes. This is simply the sum of the 
thermal entropies of the individual modes of every field specie 


a=N n 


a =1 k a =1 


k =1 


-/3pfc 


S = £ £ UK) = iv£ U T—5T - lo e(i - ^ Tk ) 


(3.56) 


where we remind that /3 = Due to the simplicity of the previous relation, the 

critical momentum p n for which the thermal entropy of the reduced subsystem is precisely 
half of the total entropy does not change when considering N fields. Hence, the critical 
momentum dividing the QFT in two halfs maximally entangled with each other is again 


p CTit = T log 2, 


(3.57) 


giving a different interpretation to temperature in the QFT. 


4 Conclusions 

The use of random unitary dynamics is having a strong impact in the physics of black 
holes and in the context of quantum thermalization, see the nice recent reviews mm 
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and references therein. Probably the most famous example comes when considering finite 
dimensional Hilbert spaces, and when using the Haar ensemble of random unitary ma¬ 
trices. This example was first studied in [13], in which the typical properties of reduced 
density matrices were derived. Later, in [5], their results were used to compute the average 
entanglement entropy of a given subsystem. 

In this article we have taken some steps towards generalizing the previous random 
unitary framework, so that it can be applied to any physical system and phase of itf^J We 
formulated the random unitary framework in Fock space. If unitary evolution constrains 
the state of the system to live in a specific subspace of the Hilbert space, defined by some set 
of charges, the typical properties of the subspace can be studied using the simple Gaussian 
relations (2.2) and (2.3). This Gaussian ensemble is not precisely the Haar ensemble, but 
deviations between both of them are exponentially suppressed. Indeed, just with those 
Gaussian relations we were able to rederive the old results of in Appendix [Aj 

Formulating the problem in Fock space makes the study of reduced dynamics very 
natural. Genetically one can choose any subset of particles, charcaterized by several quan¬ 
tum numbers, and integrate out all the others. The main results of the article are equa¬ 
tions (2.12) and (2.13). These are generic formulas for the typical reduced density matrix 
and the typical entanglement entropy associated to it. They are valid for any desired 
subsystem in any theory. They are written just in terms of ‘constrained multiplicities’ 


--(At, Q >p, i a , i;,, 


i.e. 


the number of configurations with a given set of constraints. 


In this way we connect modern approaches to quantum thermalization concerning quan¬ 
tum information theory, and in particular quantum entanglement, with more traditional 
microstate counting apporaches. 

The standard deviations of both formulas are computed in Appendix [Aj For the case of 
entanglement entropy we developed a new method to compute the deviations. The method 
directly connects the differences between Page’s curve and the exact thermal answer to the 
well-known Wigner’s semicircle law for the eigenvalue statistics of the ensemble of Gaussian 
Unitary Matrices (GUEf 4 } see Appendix [Aj This is very satisfactory, since it connects in 
a transparent way the more traditional approach to quantum chaos provided by random 
matrices and gaussian ensembles, see m for a nice introduction to the subject, to more 
modern approaches based on entanglement entropy. 

Finally we applied the generic formulas for a massless scalar held theory in a finite 
segment, and for its generalization to the case of N held species. In this context we arrived 
at the following results: 


13 Similar approaches have been considered in |14f mmmm in the context of black holes. In 
particular [15j also considers the application of Haar ensembles to microcanonical subspaces. 

14 See [Ml for a nice introduction to the field of random matrices. 







• We showed the emergence of Gibbs dynamics with the correct effective temperature as 
a function of the total energy of the state. Any desired temperature can be produced, 
just by varying the total energy. The equation of state of the QFT is seen as the 
typical macroscopic configuration of the underlying microscopic quantum state. 


• The typical density matrix is not equal to the thermal density matrix, but to the 
microcanonical density matrix. Both are seen to be equal just in the strict thermody¬ 
namic limit. This is a huge difference with the old model described in |13j . in which 
typicality is equal to thermality. 

• The deviations from thermality are of 0{1/S) instead of 0(e~ s ). These are much 
bigger than expected, and it seems just due to energy conservation. This result 
might imply that information from black hole evaporation can be extracted within 
perturbation theory. 


• We find exponentially suppressed entanglement entropies in the random state. The 
smallest ones are of (D(Se~ s ). In the context of holography, we do not expect this to 
be captured by geometrical bulk quantities, such as geodesic lengths, since these can 
only measure entanglement entropies with a minimum size of 0( 1). 


• Drawing the analogue of Page’s curve in the QFT scenario, as a function of the 
energy scale, provides a sort of entanglement running for the random state. The 
critical momentum dividing the QFT in two halfs maximally entangled with each 
other is given by p cr i t = Flog2, giving a different interpretation of temperature in 
the quantum theory. 


We want to finish with some outlook. Since the generic equations (2.12) and(2.13) 
apply to any theory with a Fock basis structure, such as a QFT or certain spin systems, 
there are two interesting avenues for future research in this regard. 

The first one is to apply the formulas to integrable systems, and see whether Generalized 
Gibbs Ensembles can be reproduced naturally, and what are the deviations from them. This 
might ultimately provide some more insight into generic differences between integrable and 
non-integrable quantum theories. Besides, although Generalized Gibbs Ensembles seem to 
fail in some situations, see [7], the random unitary ensembles we used in this article may 
represent better the physics of the integrable system. 

The second one is to apply the framework to holographic held theories, beyond what 
is done in Section 3.2. The previous formulas seem a rigorous framework to study en¬ 
tanglement between infrared and ultraviolet domains at Enite temperature. An exciting 
possibility is that it could turn into a workable route to extract physics from the near 
horizon region of quantum black holes. Aspects of these near horizon regions might be 
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encoded in the structure of the Fock space entanglement of the field theory. At any rate 
we expect this framework to give more insights into the connections between entanglement 
and quantum gravity 01201- 
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A Average entanglement entropy vs entanglement en¬ 
tropy of the average 

One of the main objectives of this article has been to give a generic expression for the aver¬ 
age entanglement entropy of a given reduced subsystem in a random state. Mathematically 
we needed to compute the following average 

^typical = ^ A j = J dpS A = _ j dpTr(p A log p A ) . (A.l) 

We will now develop a novel method to derive the average. We will show that the method 
is able to reproduce Page’s result as a specific case, but it is otherwise generalizable to our 
QFT framework as well. 

Let us consider first the classical analogue of the previous problem. This will not be 
just a pedagogical exercise, since we will need the results of this classical case for the more 
complicated quantum one. 

Consider a random probability distribution pi, for i — 1, • * • , n. With this we mean that 
each entry p t is an independent random variable with certain probability distribution itself, 
with the only constraint of global probability normalization YliPi = 1- The average over 
the randomness in the probability distribution itself will be termed by [•••]. For example, 
the moments of pt are 

\p?} , (A-2) 

for m = 1, 2, • ■ ■. The mean will be denoted simply by \j),] = pi. The fact that they are 
independent random variables implies 

\PiPj\ = \Pi}Sij • (A-3) 

Given that the mean of a probability distribution cannot be zero, it is more interesting to 
work with the random variable 6pi = Pi~ pi, the moments of which can be easily obtained 
from the moments of p t . 

In this context, problems naturally appear when we want functions of p*. The average 
value of a given function /(p;), with respect to a single realization of pi, will be represented 
as [f(Pi)\ P 

i—n 

[f(Pi)]p = '^2Pif(Pi) • (A.4) 

i= 1 

But since the probability distribution itself is a random variable, we have that [f(pi)] p 
is also a random variable, and we are more interested in its moments rather than in the 
variable itself. In particular the mean is given by 

i=n i=n 

[[f(Pi)]p\ = E Piffa)] = Efe + Pi)f( S Pi + Pi)l ■ (A-5) 

2=1 1=1 
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To compute the average of the function we can proceed iteratively, by Taylor expanding the 
function in terms of the deviations from the mean. For example, in the case f(pi) = — log Pi 
we are computing the average of Shannon’s entropy: 


[/(Pi)]p = S(p) = -J^pilogpi , 


(A.6) 


which can be Taylor expanded as 


n 


S(p) = ~^2^Pi+Pi)^g(5pi+pi) 


1=1 




(A.7) 


Since [dpi] = 0 by construction, up to higher order terms in the Taylor expansion, the 
average Shannon’s entropy is given by 



(A.8) 


where DEA stands for Deviation from the Entropy of the Average, a quantity that needs 


to be computed to make sure the first term is the leading term in the themodynamic limit, 


and for other potential applications as well. 

We want to remark here, as it is more transparent now, that DEA is not a priori related 
to deviations from some definition of thermality. We want to make the case that there is 
a conceptual, precise and computable, difference between DEA, the Standard Deviations 
from Typicality a 2 (S(p)), and Deviations from Thermality DT. Indeed one could also be 
interested in the average deviation of a typical realization from the mean [5'(p)|^| This is 
the standard deviation from typicality 


° 2 (S(p)) = [£(p) 2 ] - [S(p)] 2 . 


(A.9) 


To compute it, we need to find the first term in the sum using again the previous expansion 
method. We obtain 



(A.10) 


15 This quantity has not been studied before in the context of quantum thermalization. We comment on 


this below. 
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The average of this expression is given by 


[5(p) 2 ] ~ S(p) 2 + ^[(Spi) 2 } f (log pi) 2 - 


Sjp) 

Pi 


(A.ll) 


and combining formulas we finally obtain: 

& 2 {S(p)) = ^[(^i) 2 ]( lo gPi) 2 + ■ ■ ■ , (A.3.2) 

i 

where the dots indicate higher order terms in 5p. 

We have written the average entropy (A.8) and the standard deviation from it (A.12) 
in terms of \{8pi) 2 ]. If we could obtain this quantity for the problem of interest we could 
obtain the average quantities. 

Let us come back now to the quantum case. As with the classical case we will call 
[p] = p. For all cases discussed in this article we have (p)ij = (p*)ij = ( p)u8ij , i.e the 
typical density matrix is a real and diagonal matrix. We follow the same procedure and 
change the integration variable, from the state p to the deviation from the mean: 


p = 5p + p , 


so by construction we have 


[(<*p)«'] = 0 > 

where %%' are the matrix indices, and also 

[(<5p )ii'(fip)jj'] [(P — (p) ii (p) jj^iV^jj' 


(A.13) 

(A.14) 

(A.15) 


We conclude that the reduced density matrix will be some typical density matrix plus some 
random matrix, with a structure that might depend on the case studied. 

Finally, we can translate the average over p to the average over its random eigenvalues 
Pi- 

[5 a ] = ~ f dpTi(p A log p A ) = [-^Pilogpi] , (A.16) 

i 

and therefore, at first order, over the random eigenvalues 5pi of the random matrix 5p: 

[5 a ] = [- y1(pi + Spi) log (pi + 5pi)] . (A.17) 

i 

In this way, if we know the statistics of the random matrix 5p and the statistics of its 
eigenvalues, we can compute the average value of the entanglement entropy and its standard 
deviation using formulas (A.8) and (A.12). We will compute it for bipartite systems next. 
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A.l Random dynamics, GUE ensembles and Page’s formula 


Following Ref [5, 13] we consider a bipartite system, with total dimension AB, where A is 
the dimension of subsystem A, and B of subsystem B. If |z), with i — 1, • • • , A, is a basis 
of A and |a), with a = 1, • • • , B, is a basis of B, the state of the system can generically 
be represented as 

10 =£*.!<,«> (A.18) 

i,a 

The application of a random unitary to the previous state can be defined as in Section 2, 
except that we let it act in the full Hilbert space instead of the subspace of fixed energy. 
It reads 

[V’ia] = 0 bl’jp'l’ia] = ASijSap . (A.19) 

At first order, the random unitary just produces AB independent gaussian random vari¬ 
ables with squared deviation given by A. To fix this deviation we just fix the average 
normalization of the state 


(WM = y [CA«] =A = l. (A.20) 

i,a 

We can now trace out subsystem B. The reduced state for subsystem A is given by 

P A = • (A-21) 

i,j a 


We can now derive the statistics of the reduced state from the previous statistics (A.19). 
The average is given by 

(A. 22 ) 


\P A \ = = 4 A > 


so the average density matrix is exactly equal to the thermal density matrix at infinite 
temperature. The mean probabilities are given by pi = 1/A for each i. We remark here 
that this is not generic, as we have shown in the article. When charges are conserved the 
average density matrix is naturally the microcanonical density matrix, and therefore it is 
close to the thermal density matrix, but not equal. In such cases computing the deviations 
from the entropy of the average is not as important, because the leading behavior of the 
deviations from thermality are just given by using the entropy of the average. 

The deviation from the mean is seen to be 


\{P A ) 


r n A 

Pjj' 


a,(3 






Ay + 


A* B 5ij5i ' j ' 


(A.23) 


Therefore, the statistical properties of 5p A = 
more significantly 




p — p A are just given by [(£p A )r,-] 
= ~A^B^ ik ^ jl ' 


0 , and 
(A.24) 
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We conclude that the statistics of 5p A are those of a random matrix belonging to the 
Gaussian unitary ensemble (GUE), with o 2 (5p A ) = Given this insight we can use 

one of the basic results in random matrix theory, see [23=] for example, which is the single 
eigenvalue probability distribution /(A). The function /(A) is the celebrated Wigner’s 
semicircle law. For a matrix of size N , with deviation a for each of the entries, it is given 
by 

/(A) = 


V4AW 2 


A 2 , 


(A.25) 


47riV<7 2 

so that [A] = 0 and [A 2 ] = No 2 . In our case N = A, and this implies that the eigenvalues 
of 5p A , given by 5pi, satisfy 


\m 2 \ = 


AB 


Using now relation (A.8) we obtain: 


[Sip)} = ~ J2 Pi logPi ~ l 

■ ^ ■ -i Pi 

i=l i=l 


2 ) = s(p)-A, 


(A.26) 


(A.27) 


which is, to the level of accuracy in which we are working, Page’s result for the average 
entanglement entropy [5j. The deviations from the entropy of the average are then 


DEA = — . 
2 B 


(A.28) 


When A ~ Oil), the DEA is inversely proportional to the size of B, or equivalently, 
exponentially suppressed in the entropy S = ln(AB). On the other hand, when A ~ B, 
the DEA is of C>(1). 

In this case, because the mean density matrix is exactly equal to the thermal density 
matrix, these DEA deviations are equal to the deviations from thermality DEA = DT = ^A, 
a feature which will not extend to the QFT case. 

A nice feature of this framework is that it can be immediately extended to compute 
the deviations from typicality, which are given by relation (A. 12) 


O' 


; (S(rt) = ^[(fe) 2 ](iogB) 2 + ■ • • = 


(A.29) 


The standard deviation from typicality is therefore parametrically smaller than the previous 


DEA =4 


A.2 Average entanglement in a QFT framework 

The extension of the previous framework to the QFT case does not entail any more con¬ 
ceptual insights. It just brings some technical complexity, because the structure of the 
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reduced density matrix is not as simple as the previous case. But the final result turns 
out to have the same form, i.e we will obtain DEA = A, for an appropriate A and B that 
we determine below. What is different is that for the QFT case we have DEA ^ DT, 
since the typical reduced density matrix is not equal to the thermal one. Therefore we 
rigorously justify the use of the entropy of the average through the article, since it provides 
the leading difference from thermality in scenarios in which several charges are conserved. 

Let us first consider the case of a reduced density matrix for one particle type a. 
This case is particularly easy since the reduced density matrix is already diagonal before 


averaging, see (2.7). The diagonal entries, or the probabilities, are given by 


Via = 


'la,OL 


(A.30) 


We can therefore compute the DEA directly from (A.8). One easily finds 

Q(E t , i a ) 


Via = 


[fej 2 ] = 


n(Er) ’ 

The deviation from the entropy of the average is then 


0 (E t , i a ) 

0(£ t ) 2 


1 (V m 
DEA a = - ' ° 

2 


+ 1 ) 


2 * 


(C ax + l)e 


-s 


(A.31) 


(A.32) 


0 (£ t ) 

hence exponentially suppressed in the entropy S. 

Now we generalize to any desired subset of particles a,b, ■ ■ ■ , c. This is more difficult, 
since the reduced density matrix is no longer diagonal before averaging. To derive the 
deviations from the entropy of the average we have to understand the structure of the 


reduced density matrix. The reduced density matrix is given by formula (2.11), which we 
rewrite here 


Pa,b,-,c = ^<«|'L ou t)(T out |a) (A.33) 


a 



where a labels the set of particles which are traced out. Although this is not a diagonal 
density matrix we notice that it has a certain block diagonal structure, due to energy and 
charge conservation. The entries outside the blocks are direcly zero, without the need of 
averaging. The blocks are defined by their total energy E^’ lb) '" ,lc and their total charge 
Qia,i xhe dimension of each block is, by definition, {l(E^’ lb ’"' ,tc ,Ql£’' lb ’"' ,%c ), another 
multiplicity that can in principle be computecp*} 


16 The generating function of such a multiplicity is indeed easy to write with the techniques described in 
Appendix [Bj But we will not need its exact expression in a given theory, as we show below. 
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The statistics of the reduced density matrix are its average 


Pa,b ,c — [ Pa,b ,c] ^ ^ 




Qt) 

and the deviation from the mean, whis is given by: 


|fa; ffc; i ^cji^ai ib] 


(A.34) 


[(p)*Ap)jj'} = 


Q(E T ,Q^,if 


Oii'Oii' T 


fi(Er, Q r T , i ) 


On Or 


n(E T ,Q r T y n n n(E r ,Qi p ) 2 y ’ 


(A.35) 


where we have used a short hand notation for i = - ■ ■ ,i c , etc, and where it is 

implicitly assumed that i,i!,j,j' belong to the same block, since in a different case the 
entries are directly zero. 

From the previous expressions we can compute the statistics of Sp a ^... jC = p a ,b, -,c ~ 
p a ,b,-,c■ By construction [8p a ^... )C ] = 0. The deviations are given by: 


l(fipa,b,— ,c)ii' (8pa,b,— ,c)jj'} ~ 


r2(A r r, Qti i*at ^bi j g) 
tt(E T ,Q r T ) 2 


OijOj/'j 1 


(A.36) 


We conclude that 5p a} b,-, c is a block diagonal matrix, with each block Bi characterized 


by its total energy E, 


1 a fib j" ” i^c 


and its total charge Q\ 




, and where each block is a 


random matrix taken from GUE with size Sb, = ,lc , ,ic ) and deviation 


2 _ " j^c) 

a Bi ~ n(Er,Q^) 2 ' 


If Spi are the eigenvalues of the block Bi, the contribution to the deviation from the 


entanglement entropy of the average of Bi is (see formula (A.8)) 


DEA ®. = 5 E 

Pi 


Sl 


m£Bi 


2 Q{Er,Q r T ) 


(A.37) 


The total deviation between the average entropy and the entropy of the average is finally 

(A.38) 


dea = •_ 

2 Qj 


Although this expression seems opaque at first sight, it is exactly the same expression as 
for the Page’s case. Notice that JA S% counts the number of non-zero terms in the reduced 
density matrix, which for one particle subsystems is just z™ ax + 1. There is an analogous 
expression for the reduced density matrix of the complementary subsystem, say ^A S 2 R . 
The direct product of these two numbers must be the total non-zero entries of the full 
density matrix, so we have JA S%. JA S 2 B , = £1(Et, Qt) 2 - We conclude that y^JA S^. = A 
provides the right dimension of the c hosen subsystem, while for the complement ary w e 
have \/Y2i^B- = B. Relation (A.38) is therefore exactly equivalent to relation (A.28). 


Finally, the deviations from typicality can be computed using (A.12) in a similar fashion. 
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We conclude that the difference between the average entropy and the entropy of the 
average density matrix behaves in the same way as the common Page’s case [5], i.e it is 
of 0(e~ s ) for small subsystems with dimensions of 0(1), and it is of 0(1) for subsystems 
with dimensions of O(S). 

On the other hand, the deviations from thermality are naturally defined by the differ¬ 
ence between the thermal entropy and the average entanglement entropy 

DT = S 0 - [5s] ~Sf,-Sp + DEA . (A.39) 

In Page’s case, the average reduced density matrix is equal to the thermal density matrix, 
and therefore the first two terms cancel each other, leaving DT = DEA = A/2B. In the 
QFT case, since for small subsystems we already have Sp — Sp ~ (1/S), we do not need 
to consider the subtle DEA, since the leading corrections are already accounted by the 
entropy of the average density matrix. 


B Partitions and constrained partitions 


One of the main messages developed in the article concerns the direct connection between 
typical reduced dynamics and multiplicity counting. In particular the computation of 
reduced density matrices and their associated entanglement entropies boils down to the 
computation of microcanonical and constrained degeneracies. In the specific cases consid¬ 
ered in this article this turns out to be possible, and we provide the detailed process in 
this appendix. 

In the case of the real scalar field on a lign segment, we need to compute Q(Et), the 
number of states at energy E- p, and VL(E^,i n ), the number of states at energy E-r with i 
particles with momentum p n . f2(-&r) was seen to be the number of different partitions p(n) 
of the natural number - E t, which is given asymptotically by 


(l(E T ) = p(-B T ) -> _JL_eVW . (B.l) 

7T 4V3 LE t 

In the same way, fl(E^,i n ) is the number of different partitions p(n',i n ), of the natural 
number n' = ^ E t, in which n appears i n times. To compute it, notice the following 
theorem due to Euler: 


P(x) 


111 
1 — xl — X 2 1 — X 3 


1 

1 — X 1 


Y p(m)x m , 

m =0 


(B-2) 


which provides a generating function for the number of different partitions p(m). This is 
because 


P(x) — (1+x+x 2 -] - )(l+a; 2 -l-a; 4 -|-)(l+a; 3 -l-a; 6 -| -) ■ ■ • (l+x n +x 2n -\ - )■•• , (B.3) 
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and picking a monomial in the k-th part looks like x lhk for arbitrary positive integers %k- In 
the product P{x), we get all monomials of the form x lll x t22 x 133 ■ ■ ■ = x n+2 * 2+3 * 3 +"' , which 
produce all possible partitions of the integer 

m — ii + 2 i 2 + 3i 3 H- ni n , (B.4) 


with n < m and i m — 0 or 1 . 

To fix the number of times a given number n appears in the partition, we just need 
to fix the monomial in the corresponding parenthesis, i.e we fix and select in the n-th 
bracket only the monomial x lnn . This way, we get the partition of m in which the number 
n appears i n times for fixed (n, i n ). Hence, the generation function of p(m, i n ) is given bjQ 

Qiu( x ) = (1 + x + x 2 H- )(1 + x 2 + x 4 H- )(1 + x 3 + x 6 H-) ■ ■ ■ 

x (1 + x 11 - 1 + x 2{n ~ l) + ■■■ )x inn (l + x n+1 + 

m= oo 

= p(m, i n )x m . (B.5) 

m =0 

In fact, the sum starts only from m = i n onwards, so p(m, i n ) = 0 for m < i n , which is also 
obvious from the definition of the partition. The previous expression implies 

m— oo 

Q in (x) = x inn (l - x n )P(x) = {p(m)x m+inn - p(m)x m+n ^ +1) ) , (B. 6 ) 

m=0 

so we get, as long as m — i n n > n, 

p(m, i n ) = p(m — i n n) — p{m — i n n — n) . (B.7) 


For m — i n n < n, the identity is still true but without the second term. We now finally 
obtain 

Q(E T ,i n ) = p(-E T ,i n ) = p (-E T - i n n) - p (-E T - n(i n +1)) , (B.8) 


7r 


n 


71 


a result that was used in Section 3.1. 

One of our main result in the text is the sub-leading corrections to the thermal entropy, 
given in (3.20). It is easy to proof this formula using the Hardy-Ramanujan asymptotic 
formula for the number of partitions, 

P(N ) 


11 ^ PIN 
-e v 3 


t ,/3 A' 


(B.9) 


17 The partition p(m,i n ) should not be confused with the partition denoted by p(m,n), which is the 
partition of m with largest part n. 
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Consider now two numbers a and b, with a < b and both a <C N and b -C N. Define the 
combination 


P(N , a, b) = 


a ) — p(N — b) 


P(N) 


(B.10) 


which is of the form of the reduced probabilities in (3.15). Using (B.9), it is now an easy 
exercise to show that, in the a/N expansion, we get 


d^d= s (l + ^ + ...)e-^( 1+ « + "0 . (B.U) 

Higher order corrections are of order a 2 /N 2 or more. One would be inclined to also drop 
the subleading term in the exponent, such that the only next to leading order term is the 
a/N term in front of the exponent. This would be fine, unless a ~ \fN, which is what we 
actually consider in the main text (when the momentum mode has energy of order T). In 
that case, a/N ~ l/y/N is still small, but expanding the exponential now gives subleading 
terms that are of the same order as a/N. The terms don’t cancel out each other but merely 
change the coefficient in front of a/N. So we conclude that 


P(N, a, b) 




_ 7T O' 

e n 


1+ °N 


!) . \ _ b 


7T b 

e v&n 


(B-12) 


The corrections are hence of order a/N, and if we take a ~ vN (and similarly for b), 
the exponent is of order unity, and the total error scales like 1/y/N. In the main text, 
1/Vn 1/S , so the error is inversely proportional to the entropy. 

Generalizing this procedure to include more constraints is straightforward. We could 
be interested in ■ ■ ■), the number of states with i n units of momentum n, ii 

units of momentum l, etc. This is then equal to p(m,i n ,ii , • • •) the number of partitions 
in which n appears i n times, l appears i\ times, etc, and where m = We get this 

number from its corresponding generating function: 


- x n )x il \l -x 1 )--- P{x) . (B.13) 


The previous formula can be applied for example to compute the entanglement between 
high energy momentum modes and low energy momentum modes. If we trace out momenta 
bigger than p n , we need to analyze the following generating function: 

Qh,i 2 ,-,in( x ) = x ll (l - x)x l22 (l - x 2 ) ■ ■ ■ x ln1l (l -x n )P(x) . (B.14) 


To convert this to a generating function, we use the following identity: 

n 

(l-x)(l-x 2 )---(l-x n ) = l-J2 xl + J2 xl+r - x l+r+s + • • • + x 1+2+ - +n , (B.15) 

1=1 l<r l<r<s 
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where e.g. 

n— 1 n n —2 n —1 n 

E = £E. E=EEE- p-w) 

Z<r Z=1 r=l -\-1 l<r<s 1=1 r=l -\-1 s=r+l 

This provides a very complicated expression for the number of modes with q particles with 
momentum p\,i 2 particles with momentum P 2 , etc., which is given by: 


Q(E T ,i u i 2 , ■ ■ ■ ,i n ) = p ( — E t - y^j k k j + 


k =1 


- 1L,p[^ E r-J2 lkk - 1 ) + 


«=i 


k= 1 

n 


l<r \ k =1 


(B.17) 


fc=l 


(=1 


Although this seems an opaque expression, it can be matched exactly with predictions from 
the Gibbs distribution, as was explained in a previous section. 

Now we describe the computations in the second case, a theory of N scalar fields. 
We will proceed in much the same fashion as the previous case. First we compute the 
generating function Q N (x) of the number of ways p N (n), with n = LE^/n, of writing: 


a=N n=n n 


n = 


^ in a n a t 


(B-18) 


a=l n a = 1 


It is simple to observe that: 
Q N (x) = ^2p N (n) 


x = 


- 1 \N( 1 \N( 1 \N 

1 — x 1 — x 2 1 — a: 3 v 1 — x i 




\N 


(B.19) 

This is just seen by expanding each of the fractions, and verifying that the products provide 

(B.20) 


all terms in (3.3). Therefore we obtain: 


P N (n) = Y P( m )p ( r )' ’ ■ p(l)5 n , m+ r+...+i , 

ra,ry-- ,Z 


where the number of indices m,r, - ■ ■ ,1 is equal to N, and p(n) is the usual number of 
partitions. The previous sum can be evaluated by a saddle point approximation, when n 3> 
N. Physically, the saddle point approximation is just the statement that with overwhelming 
probability the total energy Et Ntt/L will be equally distributed over all N scalar fields. 
Mathematically we have m = r = ••■ = / = n/N , and the previous degeneracy reads: 

p N (n) ~ (p(n/N)) N -> \ N e 2 *^? , (B.21) 


41 







where A = 


In this n N limit, the entropy now provides the usual Cardy formula: 


S ~ 2n\I -cn 
6 


(B.22) 


Finally, the generating function Qf qa ( x ) of the number of states of total energy E -r and i qa 
units of momentum q a associated to the held a is given by 




\N 


1 — X 1 — X 2 


1 — X q ' 




\N 


1 — X q+1 ' 


so that 


Pi q S n ) = p( m )p( r ) ■ ■ ■ P( l ) S n,m+r+---+l+i qa q 0 ~ 

,1 

~ P(m)p(r)---p(l)S n,m+rH H+q a (iq a +l) ’ 

,1 

In the previous n N limit, the saddle point approximation results in: 

,n~q a (i qa + l). 


Pf qa (n) = ( P( 


,n -i qa g a . 
N 


P(- 


N 


(B.23) 


(B.24) 


(B.25) 
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